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The one-step model of valence-band photoemission and inverse photoemission from 
single-crystal surfaces is reformulated for generalized (non-local, complex and energy- 
dependent) potentials. Thereby, it becomes possible to account for self-energy correc- 
tions taken from many-body electronic-structure calculations. The original formulation 
due to Pendry and co-workers employs the KKR multiple-scattering theory for the calcu- 
lation of the initial state. This prevents a straightforward generalization of the one-step 
model to non-local potentials. We therefore consider the Dyson equation which is set up 
within a muffin-tin-orbitals representation as an alternative to obtain the initial-state 
Green function. This approach requires a revision of the transition-matrix elements 
which is carried out in detail. The final state is considered as a time-reversed LEED 
state as usual. The proposed generalization of the one-step model allows to distinguish 
between the bare photocurrent reflecting the (quasi-particle) band structure and the 
secondary effects due to the (dipole) selection rules and due to the wave-vector and 
energy dependence of the transition-matrix elements. 
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I. INTRODUCTION 



Much information on the electronic structure at crystal surfaces is gained by ultravio- 
let photoemission spectroscopy (PES) and inverse photoemission spectroscopy (IPE) 
[P-[TO|| . Energy- and angle-resolved PES and IPE are useful tools for analyzing the dispersion 
of (quasi-particle) bands. For the sometimes rather involved interpretation of measured spec- 
tra, a comparison with theoretical results for the electronic structure is inevitable. However, 
measured energy- and angle-resolved spectra can hardly be compared with the calculated 
band structure directly: Secondary effects due to the wave-vector and energy dependence 
of the transition-matrix elements and due to selection rules considerably distort the "bare" 
spectra. Even if the primary interest rests on the bare quantities, such as the density of 
states and the dispersion of bands, a theory of photoemission is needed to build the bridge 
between the photoemission experiments and the electronic-structure calculation. 

The so-called one-step model of photoemission has quite successfully fulfilled this demand 
in recent years. The original method as proposed by Pendry and co-workers |Tl]-|13| has been 
the starting point for several generalizations and improvements, which became necessary as 
the experimental technique was more and more refined and growing interest was spent on 
more complex materials. Inverse photoemission in connection with a realistic model for the 



surface barrier |14,15.q] has enlarged the applicability of the one-step model. The relativis- 



tic generalization |l6|-p0|] made possible investigations of spin-polarized photoemission from 
non-magnetic materials . The extension of the one-step model to magnetic materials 

resulted in theoretical studies of the magnetic circular and linear dichroism [2^,24]. Pho- 
toemission from covalently bonded systems, adsorbate-covered surfaces and from materials 
with a relatively open crystal structure brought the need to overcome the muffin-tin po- 
tential model on which the original one-step formulation was based. The generalization of 
photoemission theory with respect to space-filling potentials [^,0] and the development of 
a relativistic, full-potential one-step model for materials with several atoms per unit cell 
has more or less completed recent work on that matter. 

The one-step model of (inverse) photoemission starts from a given electronic potential 
which is needed for the necessary construction of multiple-scattering states and propagators. 
This potential is taken from a separate electronic-structure calculation and must be con- 
sidered as the decisive input for the one-step model: The peak positions in the calculated 
(inverse) photoemission spectrum directly refiect the energetics of the electronic structure 
that evolves from the potential given. The reliability of the PES (IPE) theory is thus inti- 
mately related to the "quality" of the potential for the system at hand. 

The one-step model has been originally formulated for a local potential V{r). A local 
potential is provided by density-functional theory (DFT) within the local-density approx- 
imation (LDA) p8|-p3|. In practice, such LDA inputs for photoemission calculations have 
been proven to yield rather satisfactory results when compared with experimental spectra 
from several materials P^-|5B|. The use of LDA potentials has thus become a well-known 



and widely used concept in the photoemission theory. 

On the other hand, when applying the theory to strongly correlated electron systems. 
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it surely becomes doubtful whether this concept is still adequate since implicitly the LDA 
eigenenergies are interpreted as to be the one-particle excitation energies of the system. It 
is well known that there are two possible sources of error connected with that interpreta- 
tion: Firstly, the LDA only provides an approximate expression for the (local) exchange- 
correlation potential. Secondly, even with the exact exchange-correlation potential at hand, 
one is left with the problem that there is no known correspondence between the Kohn-Sham 
eigenenergies and the one-particle excitation energies PB| , |57| -P^. 

For an in principle correct description of the excitation energies the non-local self-energy 
has to be considered. This, however, constitutes a many-body problem. Therefore, DFT- 
LDA calculations must be supplemented by many-body methods to arrive at a realistic 
description of the one-particle excitations in strongly correlated systems. To give an example, 
let us mention the GW approximation |^0| which is well suited for the case of insulators and 
semi-conductors but has also been applied successfully to transition metals [^0|-|4^ . Another 
approach is to consider Hubbard-type models where those Coulomb-interaction terms are 
included explicitly that are assumed to be treated insufficiently within DFT-LDA. As some 
recent studies of this type show, significant improvement upon the LDA predictions for 
transition metals is possible indeed p4|-pT[]. 

The self-consistent potential U{r, r', E) resulting from the GW approach, consists of a lo- 
cal part ^(r) which can be identified with the LDA potential and of the non-local self-energy 
E(r, r', E) which is energy-dependent and complex. A tight-binding one-particle basis is used 
by those approaches that start from (multi-band) Hubbard-type models. Consequently, the 
resulting generalized potential is given in the form U^,^' (E) where i refers to sites, and L is 
an orbital quantum number. Again, U can be decomposed into a term V^f,^' resulting from 
the local one-particle potential in the Hamiltonian, and into the self-energy T^f^F (E) result- 
ing from the interactions. Note that transformation to the real-space representation yields a 
non-local self-energy S(r, r', E), even if it is local in the tight-binding basis, 'Ef^f'' (E) ~ Sui. 

As concerns a theory of PES and IPE, we have to face the following problem: Strong elec- 
tron correlations imply the need for self-energy corrections to obtain a realistic description 
of the elementary excitations. Consequently, a generalized potential should be considered 
to be the starting point for the actual (inverse) photoemission theory. The non-locality 
of the potential, however, causes difficulties within the original formulation of the one-step 
model: The calculation [jl2| of the one-electron initial-state Green function G{r, r', E), which 
is directly related to the (bare) photocurrent as well as the subsequent evaluation of the ma- 
trix elements is based on the Korringa-Kohn-Rostocker (KKR) multiple-scattering approach 



52,53] and is thus intrinsically limited to the case of local potentials. This seems to pre- 
vent a straightforward generalization of the theory to non-local potentials and thereby an 
application to strongly correlated systems. 

The purpose of the present paper is to demonstrate that it is possible to overcome the 
mentioned restriction. Starting from Pendry's formula for the photocurrent [jl2[, we propose 
a reformulation of the one-step model such that a general (non-local, complex, energy- 
dependent) potential U (r, r', E) can be included via the Dyson equation in the calculation 
of the initial-state Green function. For the final state the usual (layer-KKR) multiple- 
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scattering theory [p^- |55| , [TT| is retained. Our approach requires a complete revision of the 
transition-matrix elements. This is carried out in detail. Compared with earlier theoretical 
approaches to photoemission [^^ -|62|, our ansatz keeps the basic structure of the one-step 



model, which is highly desirable for means of numerical evaluation |]12 
In a preceding study 



concerning PES and IPE from Ni surfaces, a first pragmatic 
way to include (non-local) self-energy corrections in the one-step model was presented. That 
approach, however, was based on some simplifying assumptions by which the emission from 
different subbands could be treated independently of each other. This has restricted the 
range of applicability to special regions in k-space. Furthermore, the method rests on some 
presuppositions for the pole structure of the self-energy which were shown to be adequate 
for the case of Ni, but may be too special in other cases. With the present work we also 
try to overcome these restrictions and aim at an improved concept that allows to include 
arbitrary self-energy corrections in the one-step model. 

The paper is organized as follows: The next section introduces Pendry's basic formula 
for the photocurrent. The discussion of the Green function for the initial state in Sec. |ITl| 
elucidates the difficulties that arise within the original formalism for a non-local potential. 
Sec. |V| prepares for the proposed alternative evaluation of Pendry's formula which is carried 
out in detail in Sec. for the initial states, in Sec. ^ for the final state, and eventually in 
Sec. |V11| where the transition-matrix elements are considered and all partial results are put 
together. Finally, Sec. |VII1| concludes our considerations. 



II. ONE-STEP MODEL OF PHOTOEMISSION 



The starting point for the one-step model of photoemission is Pendry's formula p^j39||: 



/™s cx /F(Ei)Im(e;k||a|G+ AG+ AtG^-^le/kiia) . (1) 

All relevant information on the electronic structure around the Fermi energy and on the 
one-particle excitations is included in the "low-energy" propagator Gf^, i. e. in the operator 
representation of the one-electron retarded Green function for the initial state: 

Gt = J J \ra) G+(r, r', E^) (rV| drdr' . (2) 

cr =f , J, refers to the z-component of electron spin. Ei is the initial-state energy, i. e. 
El = ef — fi — hu, where e/ denotes the one-particle energy of the outgoing photoelectron, /i 
stands for the chemical potential, and hu is the photon energy. The low-energy propagator 
Gi^ is directly related to the "bare" photocurrent and thereby represents the central physical 
quantity within the one-step model. 

Referring to an energy-, angle- and spin-resolved photoemission experiment, the state of 
the photoelectron at the detector is written as |ejk||cr), where k|| is the component of the 
wave vector parallel to the surface. By means of the advanced Green function Ggo- operator 
representation and taken at the final-state energy E2 = ej, we have |/) = G'2o-|e/k||cr) for 
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the final state in Eq. (|l|). Furthermore, A = (e/me)Ao ■ p is the dipole operator, i. e. the 
perturbation which mediates the couphng between the initial and the final state. A has 
been taken in the electric dipole approximation which is well justified in the visible and 
ultraviolet spectral range. Aq is the spatially constant vector potential inside the crystal 
which can be determined from classical macroscopic dielectric theory. 

PES and IPE are complemental spectroscopies. At zero temperature T = photoe- 
mission probes the occupied and inverse photoemission the unoccupied part of the (quasi- 
particle) band structure. For PES this fact is accounted for by the Fermi function in Eq. 
(|TI): fpiE) = l/{exp{E/kBT) + 1). Apart from the Fermi function the ratio jpes/jipe jg 
known to depend on kinematic factors only. It is given by the energies and emission angles 
of the emitted photoelectrons and photons, respectively [|1^]. We can thus concentrate on 
PES in the following; IPE may always be treated analogously without any difficulties. 

The basic equation (|l|) for the one-step model of photoemission contains some approxi- 
mations that should be discussed briefiy: First we note that the formula is eventually based 
on Fermi's golden rule from first-order time- dependent perturbation theory with respect to 
A . This implies that it yields the elastic part of the photocurrent only. So-called vertex 
renormalizations and three-particle correlation effects are neglected which means the exclu- 
sion of inelastic energy losses and corresponding quantum-mechanical interference terms (cf. 
Refs. 



The "high-energy" propagator Ggo- is understood to be calculated for a local potential, i. 
e. assuming So.(r, r', E2) = 5(r — r')So-(r, E2) for the self-energy at the final state energy E = 
E2. This allows us to consider the final state G^^|e/k||(T) as a (time- reversed) LEED state 
as in the conventional one-step model. It can thus be calculated by the standard Korringa- 
Kohn-Rostocker technique (layer-KKR) (see Sec. [VI| ). Furthermore, the interaction of 
the outgoing photoelectron with the rest system has been neglected (sudden approximation). 
This approximation is believed to hold well for not too small photon energies. Let us 
emphasize that the sudden approximation for the final state is consistent with a local self- 
energy at i^^ = E2. We assume that many-body effects in the final state G2o-|e/k||cr) can 
be sufficiently well accounted for by the inclusion of a spatially constant, complex self- 
energy correction Sj(i?2) in the LEED calculation as usual. This weakly energy-dependent 
self-energy correction has to be distinguished clearly from the (imaginary) optical potential, 
which phenomenologically takes into account inelastic corrections to the elastic photocurrent 
3^, and from the (real) inner potential, which serves as a reference energy inside the solid 



with respect to the vacuum level [^. The inner potential, the optical potential and the self- 
energy correction for the final state can formally be included in the definition of a generalized 
energy-dependent inner potential: Vo,2 = Vq{E2) = Vor(-E'2) +^Voi(-E'2)- 

While typical values for the final-state energy E2 may be several tens of eV or even 
more above the Fermi energy, the energy of the initial state Ei lies some eV below the Fermi 
energy. It is this low-energy range that our primary interest should concern since excitations 
that involve the states in the vicinity of the Fermi energy can be of particular significance for 
typical electron-correlation effects. These excitations are described by means of the Green 
function for the initial state Gf^. 
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III. INITIAL-STATE GREEN FUNCTION 



Within the conventional one-step model of Pendry and co-workers |TT|-p!3|] the initial- 
state Green function Gf^ is determined for electrons moving in an (effective) one-particle 
potential as it is provided, for example, by DFT-LDA: 

^LDA(r, a) = K(r) + VH(r) + Kc(r, a) . (3) 

Within the DFT ground-state formalism the external core potential V^, the Hartree contri- 
bution Vh as well as the exchange-correlation potential V^c are local functions. On the other 
hand, it is well known that for an in principle exact description of the one-particle excita- 
tions one has to consider the Dyson equation for the Green function This includes 
the non-local, complex and energy- dependent self-energy. The development of reasonable 
approximations for the self-energy constitutes a fairly complex many-body problem which 
falls outside the scope of this paper. Just as the LDA potential, the self-energy must be 
assumed to be a given quantity for the photoemission theory. 

Several techniques that have been developed in the past to account for electron- 
correlation effects beyond DFT-LDA yield explicit analytical expressions for the self-energy 
or at least a numerical data set resulting from a self-consistent calculation. Either the self- 
energy is given in its real-space representation, T,^{r,r' , E), like in the GW approximation, 
for example [^0|-|43[], or in a tight-binding representation, S^^'(i?), as it is frequently em- 



ployed within many-body approaches that consider a degenerate-band Hubbard model with 



realistic Coulomb-interaction parameters P4|-[5l[]. The many-body calculation should explic- 



itly take into account the presence of the solid surface: For PES and IPE the information 
depth is mainly determined by the electron inelastic mean free path in the solid. Due to 
the small attenuation length for low-energy electrons PES and IPE are sensitive to 



a few layers from the surface only. However, all hitherto existing approaches concerning 
the multi-band Hubbard model refer to an infinitely extended, periodic crystal lattice and 
benefit from simplifications due to translational symmetry. For the single-band Hubbard 
model there are recent attempts to take into account the breakdown of translational sym- 
metry at the surface [|6^-p8[]. An extension of such real-space many-body techniques to the 



multi-band case is straightforward. 

For the following let us assume that the LDA potential Vlda as well as the self-energy 
S fully account for the presence of the surface. The local LDA potential defines the LDA 
Hamiltonian: 

^LDA(r, 0-) = + VLDA(r, cr) . (4) 

The self-energy is assumed to be defined with respect to /ilda and given in real-space repre- 
sentation for the moment. Therewith, the generalized potential that is obtained by adding 
the self-energy to the LDA potential, 

U^ir, r', E) = 6{r - r') VldaIf, a) + S.(r, r', E) , (5) 
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in principle yields the correct one-particle excitation spectrum. Now the initial-state Green 
function may be obtained as the solution of the Dyson equation which in real-space repre- 
sentation reads: 

[i^^i - /^LDA(r, a)] G+(r, r', E,) - J S.(r, r", E,)Gtir", r', E,)dv" = M(r - r') . (6) 

Within the original framework of the one-step model (E = 0) the calculation proceeds by 
interpreting the Green function as a propagator summing over all scattering paths that take 
the electron from r' to r. Applying the KKR multiple-scattering theory []52|- p5| , p!l[1 solves 
the problem to evaluate the basic formula for the photocurrent ([l|) very efficiently and has 



been adopted in the original work |[T^ as well as (with the necessary modifications) in all 
subsequent generalizations of the one-step model |T^-^^3|-^. 

A key point of this technique concerns the following expression for the Green function in 
the presence of a single atom with spherically symmetric potential known from scattering 
theory |6|,|6|j7|: 



G(°)(r,r',Ei) = -z^E,J2Mr<Wir^)yiUnyLir) • (7) 



Im 



Here r<, r> stands for the lesser or greater of |r| and |r'|, ipii^) is the regular 

(irregular) solution of the radial Schrodinger equation at energy Ei, and yim(r) are the 
spherical harmonics. The advantage to use the expansion (0) in the basic equation (P is 
obvious: One avoids to solve the Dyson equation in the form (|^), and the problem separates 
with respect to the variables r and r'. The representation (|^ is not restricted to a spherically 
symmetric potential but can be generalized for potentials of arbitrary shape as has been 
shown by Butler et al. [[70[] . 

However, the proof in (the appendix of) Ref. |7y] also shows that a further generalization 



of (|^) to the case of non-local potentials seems to be impossible. Firstly, one surely has 
to consider ipiif) ^^^^ V'i^(^) the (radial) solutions of Schrodinger's equation for a non- 
local (atomic) potential. But even with the correct solutions, the expression (0) does not 
solve the problem: The reason is that due to the non-locality of the self-energy the cases 
r < r' and r > r' are mixed and cannot be treated independently in Eq. (^) to show that 
[El — /iLDA ~ ^(-E'l)] G^'^^Ei) = for r ^ r'. On the other hand, however, it is necessary to 
work with r< and r> in (|^) in order to get the delta function on the right-hand side of (^) (cf. 
Ref. []7D[). Furthermore, one has to consider multiple-scattering corrections to Eq. (|^ when 
embedding the atom in the crystal lattice. However, there is no generalization of the KKR 
multiple-scattering formalism to our knowledge that works for a non-local lattice-periodic 
potential. 

Within the original one-step formalism the representation (|^) is employed from the very 
beginning. Consequently, for the case of non-local potentials we have to disregard (0) as well 
as the multiple-scattering formalism for the calculation of the initial-state Green function 
Gf^ in Pendry's formula. 
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IV. ALTERNATIVE EVALUATION OF PENDRY'S FORMULA 



The alternative is to solve the Dyson equation (^) directly. This, however, requires a 
revision of the formalism to a large extent. Only for the calculation of the final state the 
multiple-scattering approach can be retained. In the following it is shown that a consequent 
reformulation in the sense mentioned is possible indeed. 

The Dyson equation is most conveniently solved within a tight-binding representation. 
For this purpose we choose a set of one-particle wave functions, 

$,i,(r - R,) = {r\iLa) , (8) 

which span the subspace Tis that includes all s-, p-, d- (and f-)like eigenstates of the LDA 
Hamiltonian /ilda in the vicinity (several eV) of the Fermi energy (minimal basis set): 

E <^.Ur - R.) {S;')':f $*av(r' - R.') = 6{r - r') . (9) 

ii'LL' 

The wave functions $iLcr(i" — Ri) are assumed to be well localized at the respective sites Rj 
of the crystal lattice. L = {l,m) is a composite index characterizing the electron's angular 
momentum, and S denotes the overlap matrix: 

S^^ = I <f:L.(r - R.)<fi'L'.(r - R,)dr . (10) 

Since the basis functions are chosen to have definite angular-momentum character, the basis 
is non-orthogonal. Furthermore, the choice of the basis reflects the relevant energy range 
around the Fermi energy we are interested in. To justify the usual restriction to the subspace 
Hs, the states \iLa) are assumed to be orthogonal to all core states at all sites. 
Using the tight-binding representation for the LDA Hamiltonian, 

T;^;^' = {^La\h^r,A\^'L'a) (11) 

and for the self-energy, 

j:fJiiE,) = {tLa\J:^iE,)\t'L'a), (12) 
the Dyson equation for the initial-state Green function reads: 

E {EisfiJr: - T^: - ^mE^}G^!i'''iE^ = ^^^hu , (is) 

i"L" 

Within the tight-binding formulation it can be solved numerically by Fourier transformation 
to k or k|| space and subsequent matrix inversion. More important, direct contact is made 
with numerous many-body approaches that likewise refer to a one-particle basis of localized 
orbitals and yield the self-energy in the form of Eq. ([T^). 

For the evaluation of Pendry's formula we consider the system to be built up from layers 
parallel to the surface. Because of the damping of the final-state wave field due to the 
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imaginary part of the inner potential Voi(i?2), it is sufficient to restrict oneself to a finite 
number of layers, i. e. to a slab geometry. Perfect translational symmetry is assumed within 
every layer. We label the layers by an index For each we define Rj^ to be the position 
of a reference atom which can be thought of as a local origin for the layer i^. Atoms within 
the layer i± are labeled by an index iy. Their positions are given by Rj = Rj^ + Hij_, where 
Riii denotes a vector of the two-dimensional lattice. (For simplicity we consider a system 
with one atom per layer unit cell.) 

Two-dimensional translational symmetry is made explicit in the notation of the basis 
orbitals: 

\iLa) = \i\\i±La) , (14) 
from which we construct two-dimensional Bloch sums in the following way: 

|q||i^La) = J2 e""'"^^" \if±La) . (15) 



Here qy is a vector of the first Brillouin zone, and N\\ denotes the number of atoms within 
a layer {N\\ h-^ oo). By means of Fourier transformation to the Bloch-sum basis, the Green 
function can be rewritten as: 

<^£v(^i) = ]^ E e--"^^'i|-^V)GW-'(£;0 . (16) 

II 

The operator representation of the Green function now reads: 

Gta = E hf±La) G(;]^f;v(^i) {ci\\i±'L'a\ . (17) 

The operator Gf^ is independent from the choice of the one-particle basis. Inserting into 
Pendry's formula, we get: 

I^^^^ME,)lm Y: E^«xL.(6/,k||,q||)G(;]^f;,,(E0i^il'L'.(e/ (18) 

i±i±'LL' q|| 

where M is the matrix element of the dipol operator between the final state |/) = G2a\^f^\\<j) 
and the Bloch sum \(i\\i±La): 

Mi^Laief, k||, qii) = (e/k||(7|G'+ A|q||ixL(7) . (19) 

The final state G2CT|e/k||cr) is an eigenstatc of /ilda with eigenenergy E2 = e/. Let us 
also consider the low-energy eigenstates of /ilda within the subspace Hs which we refer to 
in the following as the initial states. These may be characterized by the spin projection a 
and the parallel wave vector qy which are good quantum numbers: 

/iLDA|riq||Cr) = e„(,(q||)|nq||(T) . (20) 
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Here e„o-(q||) are the LDA eigenenergies. The LDA eigenstates form an orthonormal basis 
of Tis- We make use of them to rewrite the transition-matrix element: 



Mj^i^(ej,k||,q||) = ^(e/k||(T|G'J^A|nq||Cj)(nq||CT|q||U^or) . 

n 



(21) 



Thereby the matrix element of the dipole operator is given between eigenstates of /^lda- 
This is shown below (Sec. [VI I|) to facilitate the calculation to a great extent. 

We now proceed as follows: In the next section we will deal with the initial states 
|rzq||(T), Sec. [V^ focusses on the final state |/) until finally we are in the position to derive a 



computationally feasible expression for the transition-matrix elements in Sec. VIJ 



V. THE INITIAL STATES 

Let us "switch off" the non-local self-energy for a moment, i. e. take S = 0. In this case 
the Green function becomes diagonal in the orthonormal basis of the initial states: 

Gt = Y.\n^\\^)GU,^m{nci\\a\. (22) 

nqy 

We then have 

- hmG^^^^^Er) = M{Er - (6„.(q||) - /i)) (23) 

as a factor within Pendry's formula (0). For the case S = we can thus conclude that it 
is sufficient to consider merely a single initial state, namely the one with the eigenenergy 
e„o-(k||) = El + fi = e/o-(k||) — hij. If the self-energy is switched on again, its imaginary part 
will cause a broadening of the delta-function indicating a finite quasi-particle lifetime. In 
a sense the LDA band structure will be smeared out to some degree. Furthermore, strong 
correlations may lead to satellites in the excitation spectrum which take some spectral 
weight from the main bands. Both effects imply that for a given photon energy hu there 
may be different initial states |nq||(T) coupled to the same final state |/). Opposed to the 
conventional evaluation of the one-step model, we thus need all eigenstates of /^lda around 
the energy e/a-(k||) — hu. 

We expand the initial states in the Bloch-sum basis: 

l^qy^) = Oi'^^Lh\i±L(^) ■ (24) 

The expansion coefficients as well as the LDA eigenenergies have to be determined 

from the following non-orthogonal eigenvalue problem: 

E Klv(qii) - enM\\)St'i'M\\)) = , (25) 

which is easily derived from (|20|) when expressing the LDA-Hamiltonian and the overlap 
within the Bloch-sum basis: 
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LV 

Stt'M) = (qilU^o^lqilu'I^V) . (26) 



^ixixv(qil) = {(l\\i±La\hLBA\ci\\i± L a) 

fLL' 



The q|| -dependent (LDA-) Hamilton and overlap matrices are connected with the Hamilton 
and overlap matrices from Eqs. (p!oD and (|ll|) via two-dimensional Fourier transformation: 

II qy 

^^v = ^Ee""^'''"~''^"^^;v(q||). (27) 

II 111 

From Eqs. (|^) and (p6D we have: 

.q|,a|q||uLa) = KX)* ^f;i.(q||) , (28) 



which gives us the second factor in Eq. (^) for the transition-matrix elements. 

For the first factor, i. e. for the actual matrix element, we need the real-space represen- 
tation of the initial state InqyO"): 

vI/S;Ur) = (rinqiia) = a^-^l^TUA^ - R, J . (29) 

Here we have introduced the real-space representation of the Bloch sums [cf. Eqs. (^ and 
(El)]: 

*-_fL.(r-R.J = (r|q||.^La) = -L e^'^il^'ii $.,.(r - R,) . (30) 

» II 711 



To proceed further, it becomes necessary to specify the basis orbitals $iLo-(i") and the 
form of the LDA potential. We assume that Vlda(i') can be sufficiently well approximated 
by a muffin-tin potential, which is spherically symmetric within non-overlapping muffin-tin 
spheres centered at the lattice sites Rj and constant inbetween. The muffin-tin approxima- 
tion will greatly facilitate the calculation of the transition-matrix elements since we then 
have a simple separation in radial and angular parts. This also implies the need for a one- 
center expansion of the initial states ( pUD and thus of the Bloch sums, which in Eq. ( PP] ) are 
given in terms of the basis orbitals centered at all sites i\\ within a layer i±. For this purpose 
we may choose the basis orbitals to be (augmented) muffin-tin orbitals (MTO's) ||7T| , |7^ . Be- 
ing orthogonal to all core states at all sites of the lattice and spanning the subspace Tis, they 
have all the properties we demanded above. Furthermore, we have the desired one-center 
expansion for the Bloch sum of the MTO's of a layer z^: 

^ <;^,(r - R,J = 5.,'.,$ .L.(r - R.^O - ^ Zli;^-(q||) $,iv(r - R.,0 . (31) 



The expansion converges within the muffin-tin sphere at Rj ' = of a given layer i±' |7T 



Z is the potential-independent structure constant and $ an augmented spherical Bessel 
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function. Let us mention that the existence of the one-center expansion rests on an analogous 
expansion theorem for the spherical Neumann functions and not necessarily on the muffin- 
tin form of the LDA potential $ and $ are separated into a radial and an angular 
part: 

$iLa(r) = (l)iUr)YL{r) , 

^^La{r) = Mr)YL{r) ■ (32) 

Because of two-dimensional translational symmetry within a layer, the radial parts are 
the same for all sites i\\ within a layer and depend on the layer index i± only. A further 
specification of the radial parts is not necessary for our present purposes. 
Inserting (^) and (|32D into (|^) we get 



where we have defined 

nqii ^ nqy 

B:2^=Y.zWi^c^L'.- (34) 

Therewith we have the final one-center expansion of the initial states into spherical harmon- 
ics. 



VI. THE FINAL STATE 



For the determination of the final state 



v^l,;,(r) = (rlG^- |6^k,|a) (35) 

we consider a conventional diffraction situation with |e/k||cr) (the state in the vacuum at the 
photoelectron detector) interpreted as a source of electrons which penetrate into the crystal. 
Within the framework of the layer-KKR multiple-scattering formalism p2|- |55|JTl[] we can 
construct the whole high-energy wave field from the source |e/k||cr) and from the scattering 
properties of each layer. The result gives us the time-reversed 'axial state (^\E'[f||^(r)^ . Just 
as in the case of the initial states the main goal is to derive a one-center expansion of the 
final state into spherical harmonics within a given muffin-tin sphere at r = Rj|| + Rjj^. 

To begin with, we consider the region of constant inner potential Vo,2 between the muffin- 
tin spheres where plane waves are solutions of Schrodinger's equation. The potential Vlda 
can be assumed to be invariant under two-dimensional lattice translations r i-^ r + Rj|| . For 
any multiple-scattering event this implies that the wave vector parallel to the surface ky can 
be changed by a reciprocal lattice vector gy only. In the interstitial region the high-energy 
wave field can thus be expanded into plane waves with wave vectors 
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k± = (^-k|| + g|| , ±^2m,{E2 - Vo,2)/h' - I - k|| + g|| P ) . (36) 

We define the positive z-axis pointing into the crystal and choose r = Rj^^ to be a local 
origin for a given layer i±. Consider then a plane- wave field advancing on the layer from the 
+z side: 

^^T.K..^e^^'''~''^^' ■ (37) 



Scattering at the layer i± gives rise to a transmitted and to a reflected wave field on the 
+z and the —z side of the layer with coefficients V^^^^^ and W~^i^^, respectively. They are 
given via the transmission and reflection matrices: 

g(l 

^^s\\^±o- -"-giigji '^'^gf|j_LCT • 
g'l 

Similar relations hold for an advancing wave field from the —z side. Conventional LEED 
theory [|ll|] provides us with explicit expressions for the scattering matrices, let us cite the 
final result: 

Here we have written tiij_cr = (e^*'^"^'^ — l)/2 for abbreviation. By we denote the layer- 
and spin-dependent phase shifts which completely characterize the scattering properties of 
the spherically symmetric potential within a single muffin-tin sphere. Furthermore, A stands 



for the area of the layer unit cell, k = y2me(-E2 — Vo,2)/^ and L = (/, — m). The matrix 
Xll' = t/.,.E'4vr^('-''-'")4'H-k||R„) lV(-'R,,)(-l)(-'+-")C^x^^„ e^^n^'n , (40) 

i||L" 

which is due to Kambe [Q,^, corrects for intra-layer multiple-scattering events. Here h[^^ 
denotes the spherical Hankel function and 

Clul"=I YL{p)YL>{p)YL',(Sl)d^l (41) 

J(47r) 

are the Gaunt coefficients. The definition of the X-matrix is closely related to the structure 



constants defined via Eq. (|31|) yl 



By the matrices (pQf ) the scattering properties of all layers are known. For a given plane 
wave |e/k||(T) advancing onto the crystal from the vacuum side, they determine the whole 
high-energy wave field between all layers in the semi-infinite crystal. Its coefficients can 
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straightforwardly be found within a recursive layer- by-layer scheme [|T2|]. Since the flux 
of elastically scattered electrons is permanently reduced due to the imaginary part of the 
(generalized) inner potential, the wave field can be neglected beyond a finite depth. 

Once all coefficients of the plane-wave expansions are known, the next step will be to 
derive an expansion of the high-energy wave field into spherical harmonics centered at the 
local origin Rj^^ of a layer i±. The advancing plane- wave field at the layer i± is: 

^(0). fj.) = ^y(w+- (^""i^-^^^^ + V-- e'^'^i^'-^^A , (42) 

iV|| en \ ^ 



where VFl^,- „ and VI,,- „ are the coefficients of the plane waves on the —z and the +z side 
of the layer, respectively. From (^2]) we get the spherical wave amplitudes 



^SS = 4rE4-' (-1)" K...>r(kg) + ^g;.,.>r(k^)) (43) 



of the advancing wave field: 

+ r) = E JiMY,{n) , (44) 

L 

where ji is the spherical Bessel function. Correcting for multiple scattering within the layer 
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<lL. = T.<fLlil-^A, (45) 

L' 

we arrive at the total final-state wave field inside the muffin-tin sphere at Rj^: 

<l(r))* = ^j:^ZL.e''--^ ^.M\r -R.J)YUr . (46) 

V II ^ 

ipi^ia{r) is the radial part of the solution of the Schrodinger equation for energy E2 that is 
regular at the origin. The radial function as well as the phase shifts have to be determined 
numerically for the spherically symmetric LDA potential within the muffin-tin sphere at 
Rix- 



VII. THE TRANSITION-MATRIX ELEMENTS 

We are now in a position that allows for calculating the actual transition-matrix elements, 
i. e. the first factor in Eq. (|2T|). Again we turn to real-space representation: 

(e,k|,a|G+ Alnqp) = ^ / (^^.(r))* (AqV) vl/(;),(r)rfr . (47) 

Since [V, /ilda]- = VVLDA,cr(r), the transition-matrix element can be rewritten as: 
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(e;k||a|G+ Alnqycr) = ^ ^ , \ , / Klir))* (AoVViDA,.(r)) v&J;) (r)cir . (48) 



rUe e/^(k||) - e„^(q|| 

Here we made use of the fact that both, the final state '^^^^ and the initial states are 
eigenstates of /ilda- At this point we thus profit from the transformation of the transition- 
matrix element (p!9| ) into the expression given by Eq. (^). Since Wlba,^^^) = in the 
region of constant inner potential between the muffin-tin spheres, the integral in (|4^) now 
reduces to a sum of integrals over all muffin-tin spheres: 

VVLDA,.(r) = E 0MT(r - R,)V\/LDA,<x(r) 

i 

= J2^VUr~R^), (49) 

i 

where 0MT(r) = 1 for r < tmt and 0MT(r) = for r > tmt, and tmt is the radius of 
the muffin-tin sphere. Vicr{r) = l^o-(r) is the potential within the sphere at the site i with 
Via{r) = for r > ruT- 

Another considerable simplification of the transition-matrix elements arises from the 
two-dimensional lattice periodicity of the LDA potential. Via = Vi^a'- Since both, the final 
state and the initial states, fulfill the two-dimensional analogue of Bloch's theorem, any 
lattice translation of the form r i-^ r + Rj|| leaves the integral ( ^8]) invariant apart from 
a factor exp(z(q|| — k||)Ri||). We thus can conclude that the difference qy — ky is equal to 
a reciprocal lattice vector gy. Since we can consider ky to be fixed by the energy of the 
photoelectron at the detector and by the emission angles, the reciprocal lattice vector is 
uniquely determined by demanding qy to lie within the first surface Brillouin zone. Hence, 
translational symmetry reduces the qy-sum in Eq. (|TB|) to a single term: 

qy = ky + gy . (50) 

Referring once more to translational symmetry, we also conclude that all integrals over 
muffin-tin spheres within the same layer i±_ are equal. This leads us finally to: 



(e/k||ff|G+A|r!(k||+g||)ff)= (51) 

where the integral extends over the sphere at the local origin Rj^ of the layer i±. 

For the evaluation of the integral we separate into radial and angular parts. For this 
purpose we need the expansion of the dipole operator into spherical harmonics: 

= ^(r) - E n™(Ao)n„.(^) (52) 

and the one-center expansions of the initial states (|33|) and the final state (|46|) as derived in 
Sees. and VI. The integration over the angular parts can be performed analytically and 



yields the following angular matrix elements: 
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4vr ^ - 

DlU = ~^^0 ^lm(Ao) CimVm'li-m-m') , (53) 



where we have used the definition (|41|) of the Gaunt coefficients. The integration over the 
radial parts results in two types of radial matrix elements: 

M.^u'a = e*'"^"^ r^^^Mr)^^ir)(|>^^l'Ar)dr , (54) 

in the definition of which we included the phase shifts from Eq. (^). Since ^Ji^iaij") as well 
as and 4>i^ia{T) are regular functions at r = and since V^^o-(r) oc —Z/r for r H-i> 0, 

the integrals are well defined. 

The final formula for the photocurrent now reads: 

jP^s (X h{E,) Im Y^M.M^fM) 4S,.,'.(^i) Ml,,,^{ej,^) , (55) 

where 

/ , X ^e/i ^ 1 
Mi^i, e/,k|| = — 2^ — — ■ X 



E(":'X'')'s£'i,(kii + 



X 



\i I 'L' 



\i±' LL' J 

This completes the formalism. Combining the essentials of the previous sections, Eqs. (|55| ) 
and (^) express the result in a rather compact form. 



VIII. CONCLUDING REMARKS 

This paper has proposed a theory of photoemission (and inverse photoemission) from 
single-crystal surfaces that keeps the basic structure of the well-known one-step model but 
determines the initial-state Green function from the Dyson equation rather than by means 
of KKR multiple-scattering theory. This approach requires to reconsider the calculation of 
the transition-matrix elements which has been worked out in detail. 

In our opinion it is an important advantage of the presented reformulation that the main 
physical concepts inherent in the one-step model show up in a very transparent way: 

i) The structure of the final equation ( pSj) for the photocurrent is still reminiscent of 
first-order perturbation theory in the external electric field (Fermi's golden rule). 

ii) Consequently, one can easily distinguish between the "bare" spectrum which essen- 
tially is given by the Green function on the one hand and its modifications due to secondary 
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effects arising from the wave- vector and energy dependence of tlie transition-matrix elements 
on the other. 

iii) Due to the Fermi-function cut off (for T = 0), a non-zero intensity is found for energies 
El < only, i. e. for e/ — fio; < /x: The initial-state Green function in ( |55D describes hole 
propagation, namely propagation of the residual hole in the valence band that corresponds 
to the excited photoelectron. Since we applied the sudden approximation, photoelectron 
and hole propagate independently from each other in time. However, energy conservation 
is satisfied, and thus the photoelectron takes away the information on the hole-excitation 
spectrum. The energy dependence of the photoemission spectrum originates from the energy 
dependence of the initial-state (hole) Green function and is shifted by the photon energy 

hu. 

iv) The theory is able to account for electron-correlation effects in the initial state; in 
( p5D Gf^ is the fully interacting Green function. Final-state correlation effects (beyond LDA) 
are included only phenomenologically via the generalized inner potential Vo.2- 

v) The two-dimensional translational invariance of the surface is obvious in the notations. 
The Green function and the matrix elements are diagonal with respect to qy = ky + gy; the 
same holds for spin projection a. On the contrary, the problem does not separate with 
respect to layer and angular-momentum indices. However, the result allows to distinguish 
between the different partial contributions corresponding to hole propagation from {i^,L) 
to (u',L'). 

Apart from the advantage mentioned above, the main goal of the present study has been 
to overcome the restriction to local potentials which is inherent in the original formulation of 
the one-step model. This is a necessary precondition to be able to benefit from many-body 
theories which in general yield non-local self-energy corrections to the local LDA potential. 
Up to now these could not be used in the one-step model since they are incompatible with 
the KKR formalism for the initial state. If, on the contrary, the Dyson equation for the 
Green function in a muffin-tin-orbitals representation is chosen as the starting point, an 
unproblematic connection to the results of many-body electronic-structure calculations is 
possible. Photoemission calculations based on this new concept are intended for the future. 
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